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We shall revisit the conventional adiabatic or Markov approximation, which — contrary to the 
semiclassical case — does not preserve the positive-definite character of the corresponding density 
matrix, thus leading to highly non-physical results. To overcome this serious limitation, originally 
pointed out and partially solved by Davies and co-workers almost three decades ago, we shall 
propose an alternative more general adiabatic procedure, which (i) is physically justified under the 
same validity restrictions of the conventional Markov approach, (ii) in the semiclassical limit reduces 
to the standard Fermi's golden rule, and (iii) describes a genuine Lindblad evolution, thus providing 
a reliable/robust treatment of energy-dissipation and dephasing processes in electronic quantum 
devices. Unlike standard master-equation formulations, the dependence of our approximation on 
the specific choice of the subsystem (that include the common partial trace reduction) does not 
threaten positivity, and quantum scattering rates are well defined even in case the subsystem is 
infinitely extended/has continuous spectrum. 
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Present-day technology pushes device dimensions to- 
ward limits where the traditional semiclassical or Boltz- 
mann theory fT] can no longer be applied, and more rig- 
orous quantum-kinetic approaches are imperative 0, Q ■ 
However, in spite of the quantum-mechanical nature of 
electron and photon dynamics in the core region of typ- 
ical solid-state nanodevices — e.g., superlattices Q and 
quantum-dot structures d, 0] — the overall behavior 
of such quantum systems is often governed by a com- 
plex interplay between phase coherence and energy re- 
laxation/dephasing Q, the latter being also due to the 
presence of spatial boundaries @. Therefore, a proper 
treatment of such novel nanoscale devices requires a the- 
oretical modeling able to properly account for both co- 
herent and incoherent — i.e., phase-breaking — processes 
on the same footing. 

The wide family of so-called solid-state quantum 
devices can be schematically divided into two main 
classes: (i) a first one which comprises low-dimensional 
nanostructures whose electro-optical response may be 
safely treated within the semiclassical picture (e.g., 
quantum-cascade lasers [ill [H, [Si ) , and (ii) a second one 
grouping solid-state devices characterized by a genuine 
quantum-mechanical behavior of their electronic subsys- 
tem (e.g., solid-state quantum logic gates [13: 13) whose 
quantum evolution is only weakly disturbed by decoher- 
ence processes. 

For purely atomic and/or photonic quantum logic 
gates, decoherence phenomena are successfully described 
via adiabatic-decoupling procedures [lB| in terms of 
extremely simplified models via phenomenological pa- 
rameters; within such effective treatments, the main 
goal/requirement is to identify a suitable form of the 
Liouville superoperator, able to ensure/maintain the 
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positive-definite character of the corresponding density- 
matrix operator [l^ . This is usually accomplished by 
identifying proper Lindblad-like decoherence superoper- 
ators [171 . Il8j , expressed in terms of a few crucial system- 
environment coupling parameters [l9| . 

In contrast, solid-state devices are often characterized 
by a complex many-electron quantum evolution, result- 
ing in a non-trivial interplay between coherent dynam- 
ics and energy-relaxation/decoherence processes; it fol- 
lows that for a quantitative description of such coher- 
ence/dissipation coupling the latter need to be treated 
via fully microscopic models. 

To this aim, motivated by the power and flexibility of 
the semiclassical kinetic theory [l| in describing a large 
variety of interaction mechanisms, a quantum general- 
ization of the standard Boltzmann collision operator has 
been proposed [1|; the latter, obtained via the conven- 
tional Markov limit, describes the evolution of the re- 
duced density matrix in terms of in- and out-scattering 
superoperators. However, contrary to the semiclassical 
case, such collision superoperator does not preserve the 
positive-definite character of the density-matrix opera- 
tor. 

This serious limitation was originally pointed out by 
Spohn and co-workers [l^l three decades ago; in partic- 
ular, they clearly pointed out that the choice of the adi- 
abatic decoupling strategy is definitely not unique (see 
below), and only one among the available possibilities, 
developed by Davies in [2l|, could be shown to preserve 
positivity: it was the case of a "small" subsystem of in- 
terest interacting with a thermal environment, and se- 
lected through a partial trace reduction. Unfortunately, 
the theory was restricted to finite-dimensional subsys- 
tems only (i.e., A^-level atoms), and to the particular 
projection scheme of the partial trace. 

Inspired by the pioneering papers by Davies and co- 
workers, aim of the present Letter is to propose an alter- 
native and more general adiabatic procedure which (i) in 
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the discrete-spectrum case reduces to Davies' model [2l| . 
(ii) for diagonal states gives the well known Fermi's 
golden rule [Hjlil], and (iii) describes a genuine Lindblad 
evolution (depending on the relevant subsystem param- 
eters), even in the infinitc-dimensional/continuous spec- 
trum case, thus providing a reliable/robust treatment of 
energy-dissipation and dephasing processes in semicon- 
ductor quantum devices. 

Clearly, through our adiabatic-decoupling approach, 
different markovian approximations are generated by 
choosing different projection schemes (see below). How- 
ever we stress that, contrary to standard master-equation 
formulations (2ll.[23 , these approximations are always of 
Lindblad tvpe [isfTso that positivity is intrinsic in our 
adiabatic-decoupling strategy, and does not depend on 
the chosen subsystem. 

Moreover, our approximation scheme holds true under 
the same validity regime of the conventional Markov ap- 
proach: the so called weak-coupling limit 25], where the 
subsystem density matrix in the interaction frame moves 
slowly with respect of perturbative effects. 

In order to recall the main features of the problem, 
let us consider its general formulation based on the fully 
operatorial approach proposed in |26j. Given a generic 
physical observable A — described by the operator A — 
its quantum plus statistical average value is given by 

A = tr where p is the so-called density-matrix op- 

erator. Its time evolution is dictated by the total (system 
plus environment) Hamiltonian, that can be regarded as 
the sum of a noninteracting (system plus environment) 
contribution plus a system-environment coupling term: 
H = Ho + H'; the corresponding equation of motion for 
the density-matrix operator — also known as Liouville- 
von Neumann equation — in the interaction picture is 
given by: 



~dt 
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where denotes the interaction Hamiltonian H' written 
in units of h. 

The key idea beyond any perturbation approach is that 
the effect of the interaction Hamiltonian H' is "small" 
compared to the free evolution dictated by the noninter- 
acting Hamiltonian Ho- Following this spirit, by formally 
integrating Eq. ([T]) from to to the current time t, and in- 
serting such formal solution for p^{t) on the right-hand 
side of Eq. ([1]) , we obtain an integro-differential equation 
of the form: 



dp' 
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We stress that so far no approximation has been intro- 
duced: Equations H]) and ^ are fully equivalent, we 
have just isolated the first-order contribution from the 
full time evolution in Eq. ([T]). 

According to the standard procedure, at this point one 
should identify a subsystem of interest by, for example. 



the use of a suitable projection on the space of density 
matrices [111, [25| : although this step is crucial to study 
the weak-coupling limit, it turns out, as we shall see, that 
it needs to be considered explicitly only later on, when 
talking about positivity. For the moment, wc shall limit 
ourselves to use the projection Pq, to say that no first- 
order terms are considered in this paper, i.e., we shall as- 
sume that Po (^[H'{t), Pop]j = 0. This assumption holds 

true and is physically justified for a very general form 
of the perturbing hamiltonian [2l| , and allows us to ne- 
glect the first order term on the right hand side of ([2]), 
as it gives no contribution when projected. In any case, 
one could always assume that condition above holds true, 
by assigning the first order quantity Po{H') to the deter- 
ministic part of the hamiltonian Ho (which however could 
then be more difficult to solve). 

Hence for ease of exposition we shall now study all 
the projection-independent features of our model, safely 
neglecting first order terms. 

Let us now focus on the time integral in Eq. Here, 
the two quantities to be integrated over t' are the in- 
teraction Hamiltonian and the density-matrix op- 
erator p\ In the spirit of the perturbation approach 
previously recalled, the time variation of can be 
considered adiabatically slow compared to that of the 
Hamiltonian Ti. written in the interaction picture, i.e., 
W{t') = ul{t')nUo{t'); indeed, the latter exhibits rapid 
oscillations due to the noninteracting evolution operator 

Uo{t) — e 1^ . Therefore, in the standard (and prob- 
lematic) Markov approximation the density-matrix oper- 
ator p' in interaction frame is simply taken out of the 
time integral and evaluated at the current time t. 

Following such prescription, the second-order contribu- 
tion to the system dynamics written in the Schrodinger 
picture for the case of a time-independent interaction 
Hamiltonian Ti. comes out to be: 



dp 
'dt 



n, 
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with 



dt"H'{t') = 2 dt'Ul{t')HJJo{t') . (4) 



/C = 2 



The effective equation in ([3|) has still the double- 
commutator structure in Q but it is now local in time. 
The Markov limit recalled so far leads to significant 
modifications to the system dynamics: while the exact 
quantum-mechanical evolution in ([1]) corresponds to a 
fully reversible and isoentropic unitary transformation, 
the instantaneous double-commutator structure in ([3]) 
describes, in general, a non-reversible (i.e., non unitary) 
dynamics characterized by energy dissipation and de- 
phasing. However, since any effective Liouville superop- 
erator should describe the time evolution of p correctly 
and since the latter, by definition, needs to be trace- 
invariant and positive-definite at any time, it is imper- 
ative to determine if the Markov superoperator in ([3]) 
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fulfills this two basic requirements. As far as the first is- 
sue is concerned, in view of its commutator structure, it 
is easy to show that this effective superoperator is indeed 
trace-preserving. In contrast, as discussed extensively in 
[26i] , the latter does not ensure that for any initial condi- 
tion the density-matrix operator will be positive-definite 
at any time. This is by far the most severe limitation of 
the conventional Markov approximation. 

By denoting with {|A)}, and neglecting energy- 
renormalization contributions (26j , the eigenstates of the 
noninteracting Hamiltonian Ho , the effective equation (|3]) 
written in this basis is of the form: 

-^^^ = 2 X! [^AiA2,a'iA^Pa;a^ - 7^AiA^,a;a'iPa^,A2] +H.c. 

A'j A2 

(5) 

with generalized scattering rates given by: 
2ti 

^AiA2,A;Ai = Y^AiAi^At,Ai^(eA2 " , (6) 

ex denoting the energy corresponding to the noninteract- 
ing state I A). As discussed extensively in 26], such gener- 
alized scattering rates are obtained within the completed- 
collision limit, i.e., to — > — oo, and neglecting energy- 
renormalization contributions. 

The well-known semiclassical or Boltzmann theory Q 
can be easily derived from the quantum-transport for- 
mulation presented so far, by introducing the so-called 
diagonal or semiclassical approximation. The latter cor- 
responds to neglecting all non-diagonal density-matrix 
elements (and therefore any quantum-mechanical phase 
coherence between the generic states Ai and A2), i.e., 
P\i\2 — /ai<5aiA2j where the diagonal elements f\ de- 
scribe the semiclassical distribution function over our 
noninteracting basis states. Within such approximation 
scheme, the quantum-transport equation IjlSp reduces to 
the well-known Boltzmann equation: 

^ = 5](Paa'/a'-Pa'a/a) , (7) 

A' 



where 

Pxx' = Vxx.X'X' - ^|i^;A'l'^(eA - ey) (8) 

are the conventional semiclassical scattering rates given 
by the well-known Fermi's golden rule [23l |. 

At this point it is crucial to stress that, contrary to 
the non-diagonal density-matrix description previously 
introduced, the Markov limit combined with the semi- 
classical or diagonal approximation ensures that at any 
time t our semiclassical distribution function fx is al- 
ways positive-definite. This explains the "robustness" 
of the Boltzmann transport equation ([7]), and its ex- 
tensive application in solid-state-device modeling as well 
as in many other areas, where quantum effects play a 
very minor role. In contrast, in order to investigate gen- 
uine quantum-mechanical phenomena, the conventional 
Markov superoperator in ([3]) cannot be employed, since 
it does not preserve the positive-definite character of the 
density matrix pxiX2- 

In order to introduce our alternative formulation of 
the problem, let us go back to the integro-differential 
equation ([2]) , and let us consider the following time sym- 
metrization: given the two times t' and t, we shall intro- 
duce the "average" or "macroscopic" time T = and 
the "relative" time r = t — t'. This change of variable 
has very solid bases, as it is common and well established 
in a wide variety of contests, such as Wigner's phase- 
space formulation of quantum mechanics ,.27.] , standard 
quantum kinetics Green functions (see e.g. [1^), and 
even classical radiation theory (e.g., in the treatment of 
Bremsstrahlung) : the basic idea is that the relevant time 
characterizing/describing our effective system evolution 
is the macroscopic time T. It is now easy to rewrite the 
second-order contribution in Eq. ^ in terms of these 
new time variables: 



dT 



dT 



rc T 



W T 



T 



(9) 



In the spirit of the adiabatic approximation previously 
recalled, the density-matrix operator can be taken out 
of the time integral and evaluated at the current time 
T. As already stressed, this important approximation is 
valid when the density matrix varies slowly in interaction 
picture, that is, in the limit of weak-coupling. It is now 
convenient to replace the finite-domain time integration 
over T by introducing a corresponding Gaussian correla- 

- • - 

tion function e whose width t may be regarded as 



a safe overestimation of the so-called "correlation time" , 
which is shorter than the period of the system evolution; 
indeed, for t — to greater than the correlation time, the 
time integration may be safely extended up to infinity. To 
compare with the exact dynamics in the weak-coupling 
limit, the correlation time has to finally be brought back 
to infinity: this in turn may be accomplished by a scaling 
property of the form 

t{g)^g-^T, (10) 
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for 5 ~ 0, where g is the couphng constant, T is a fixed 
reference time, and ^ > 0. So reformulating once again, 
we assume that the density matrix moves slowly with 
respect to the correlation time t. 

Focussing on the skew-adjoint part of Eq. ([9]), i.e.. 



the so-called scattering part, which is clearly our 
main interest (the self-adjoint part is just an energy- 
renormalization term and does not threaten positivity), 
we get: 



dr 



T- -T 



,/5' (T) 



(11) 



r 



We stress how the proposed time symmetrization gives 
rise to a fully symmetric superoperator, compared to the 
strongly asymmetric Markov superoperator in [25j . 

The second crucial step in order to get a genuine Lind- 
blad superoperator for the global dynamics is to exploit 
once again the slowly-varying character of the density- 
matrix operator on the right-hand side of Eq. (|lip . 
The key idea is to perform on both sides of Eq. (fTTj) 



a so-called temporal "coarse graining", i.e., a weighted 
time average on a so-called microscopic scale, a scale over 
which the variation of p*(T) is negligible. Since in the 
small and intermediate coupling regime such time-scale 
is fully compatible with the correlation-time-scale t, we 
shall perform such time average employing once again a 
Gaussian correlation function of width 4i i-e.. 



J 



dp' 

It 



{T)=- 



2ntJ- 



dr I 



(T+r'+0,[?i'(r+T'-0,p'(T) 



(12) 



r 



Moving back to the original Schrodinger picture and com- 
bining the two Gaussian distributions, the above equa- 
tion can be rewritten in the following compact form: 



with 



dp 
dT 



C,p 



^) / dt'U\t')e 



(13) 



(14) 



looking for; indeed, the operators C are always Hermi- 
tian, and such effective dynamics is positive-definite. 

Let us finally rewrite the new Markov superoperator 
(jl3p in our noninteracting basis A, defined by the (possi- 
bly generalized) eigenvectors of Ho'- we obtain an effec- 
tive equation of motion of the form 



dpxiX2 
dt 



9 X! [^AiAa.A'jAJjPA'jA^-'PAiAJj.A'iA'iPA^Aa] +H-' 



A^ A2 



This is the genuine Lindblad-like superoperator we were with symmetrized quantum scattering rates 



(15) 



J 



V 



Ai A2 , A{ XL 



—H' , H'* 

^ AiA'j^ A2A2 



27re 



exp 



4e2 



(16) 



r 



substituting the strongly asymmetrical scattering super- 
operator given by the conventional Markov approxima- 
tion Here, h has been shown explicitly, and e = f 
is a measure of the energy uncertainty in the interac- 
tion process induced by our temporal coarse graining. 
The above scattering superoperator can be regarded as a 



generalization of the conventional Fermi's golden rule to 
the density matrix formalism; indeed, in the semiclassi- 
cal diagonal case (Ai = A2, A'^ = X'2) the above scattering 
superoperator boils down to what could be considered a 
dressed vertex-smoothed version of the Fermi's Golden 
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Rule 

(17) 

and in the limit of infinite correlation-time (e 0) the 
standard scattering rates given by the Fermi's Golden 
Rule ([5]) are readily recovered. It may be argued that 
the approximate dynamics thus found suffers from lack 
of generality, the gaussian smoothing having been put by 
hand. However, one has to look at asymptotic features 
only, as the approximation is valid in the weak-coupling 
limit e — > only: in that limit, one does not appreciate 
the gaussian smoothing, but rather the (much more uni- 
versal) asymptotic character of the correlation time PH)) . 
As such, we could say that our gaussian smoothing is 
a good representative among all the possible asymptotic 
markovian approximations to the exact dynamics, that 
guarantee a positive evolution. 

In passing, we note that the transition rates could 
be regarded as a "Quantum" version of the celebrated 
Fermi's Golden rule. This should not generate confu- 
sion: of course the Fermi's Golden Rule transition rates 
are computed according to quantum mechanical calcula- 
tions (see among thousands of textbooks), but, once 
computed, they give rise to the Boltzmann equation ([7]), 
which describes a classical Markov process '30] for clas- 
sical probabilities. Instead, the transition rates p6p do 
not describe a classical Markov process, but rather its 
quantum analog: a so called Quantum Dynamical Semi- 
group [18] for the full density matrix. 

Moreover, we would also like to note that there 
are many proposed Quantum generalization of the well 
known Fermi's Golden Rule (see e.g. [3l|), which are ro- 
bust, physically and mathematically meaningful. How- 
ever, all these generalizations consider a bipartite sys- 
tem, one of which is finally traced over, the other being 
an A^-level atom (a rather particular case). On the con- 
trary, our model is much more general and refers to a 
closed quantum system, similarly to the original idea by 
Fermi |23j . 

As pointed out previously, the whole theoretical 
scheme becomes meaningful and applicable only when a 
well-defined subsystem of interest is identified (together 
with a corresponding infinite-dimensional environment), 
so that its correlation time t can be estimated, and our 
(irreversible) semigroup dynamics can correctly describe 
the projected (but fully reversible [l3l) exact Hamilto- 
nian dynamics. As final crucial step, we shall show that 
our conclusions about positivity remain valid no matter 
how the subsystem is chosen. 

To this end, we notice that the usual partial-trace pro- 
jection, when viewed in Heisenberg picture, is of the form 

Po ■ E„ in ® S„ (^Y.n Tr{i:uBn)An^ ® 1, whcrc is 

the "environment" density matrix. An and _B„ are, re- 
spectively, generic system and environment observables. 
From this structure, it follows that the partial trace is 



completely positive (see e.g. [11]). Moreover, the pro- 
jected observables, all being of the form A (E) 1, consti- 
tute a subalgebra of the global-observable algebra. Based 
upon these two key remarks, let us now consider a generic 
projection Pq on a subalgebra X of the space of observ- 
ables in Heisenberg picture, which is also a completely 
positive map, and let {Va} be its Kraus decomposi- 
tion til], so that PqA = J2a^aAVa- Then one could 
easily verify that X is made by observables that com- 
mute with each of the Va and V^. We now observe 
that, due to its symmetry, our generator keeps the same 
form of ca. (|13p also in Heisenberg picture: by projecting 
the latter with Pq, and using the completeness relation 
J2a ^a^a = 1, One Can easily write the form for the sub- 
system's generator in Schrodinger picture, dual to the 
projected dynamics on the subalgebra X: 

a/3 a/3 

Here the "quantum transition amplitude" operators are 
given, according to ([14]), by 

Dap = VaCVp (19) 

(note that these operators are indeed ^-dependent, as can 
be seen from (ITil) ) . This Lindblad form [18'] shows indeed 
that we have obtained the generator of the completely 
positive Quantum Dynamical Semigroup we were looking 
for. Moreover, the effort we made to consider this rather 
abstract class of projection is completely justified, as, 
for example, it paves the way for a new formalism for 
Quantum Transport: suppose that P is a projection in 
our Hilbert space that identifies, say, a one-dimensional 
nano-device, and let Qi and Qr project on the left and 
right contact respectively (obviously P + Qi + Qr = 1). 
Then PqA = PAP + QiAQi + QrAQr does belong to the 
class of projections we have just studied, but the chosen 
subsystem does not come from a partial trace, nor does 
it have finite dimensions or discrete spectral properties: 
its weak-coupling dynamics needs the full power of our 
theory, in contrast with the previous ones [2ll . [2^ . 

At this point a few comments are in order. As dis- 
cussed extensively in [2^ , also for the simplest case of a 
standard two- level system — i.e., a generic quantum bit — 
the standard Markov superoperator predicts a non-trivial 
coupling between level population and polarization de- 
scribed by the so-called T3 contributions. In contrast, for 
a two- level system coupled to its environment, the pro- 
posed quantum Fermi's golden rule does not predict any 
T3 coupling term (they vanish in the infinite correlation 
time limit e ^ 0), thus providing a rigorous derivation of 
the well-known and successfully employed T1T2 dephas- 
ing model In general in fact, one could show that for 
subsystems with discrete spectra, the limit e — > reduces 
to Davies' theory [2l|. However, for subsystems with con- 
tinuous spectra, such limit is not defined, but for all finite 
collision times t > the proposed approach gives i ~ T3 , 
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so T3 contributions arc indeed present, but they become 
less and less important as the collision time I is raised, as 
it must be in the weak-coupling limit, to compare with 
the exact hamiltonian dynamics (see before). 

To summarize, we have proposed a new approach to 
modeling nowadays semiconductor nanodevices, by criti- 
cally reviewing the standard Markov procedure. Indeed, 
the latter does not preserve the positive-definite charac- 
ter of the density-matrix operator, thus leading to highly 
non-physical results. To overcome this serious limitation, 
we have identified an alternative and more general adi- 
abatic procedure which (i) is physically justified under 
the same validity restrictions of the conventional Markov 
approach, (ii) in the semiclassical limit reduces to the 
standard Fermi's golden rule, and (iii) describes a gen- 
uine Lindblad evolution, thus providing a reliable/robust 
treatment of energy-dissipation and dephasing in state- 
of-the-art quantum devices. We stress that our formu- 
lation generalizes preexisting theories significantly, as it 



gives a positive dynamics for a considerably large class 
of projections, i.e. ways to chose the subsystem, and it is 
well defined even for infinitely extended subsystems (i.e. 
with continuous spectrum). In turn, on one side this 
allows to investigate subsystems with both discrete and 
continuous spectra, a feature largely shared by mesoscale 
electronic and opto-electronic quantum devices; on the 
other side, it suggests a new way to treat electrical con- 
tacts for quantum devices, thus opening up the exciting 
possibility of a new formalism for Quantum Transport. 
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